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Major Drivers of Current Work 


• To provide computational tools as an economical option for developing 
future space transportation systems (i.e. RLV subsystems development) 

Impact on component design -► Rapid turn-around of high-fidelity analysis 
Increase durability/safety -t Accurate quantification of flow 

(Le. prediction of flow-induced vibration) 

Impact on system performance -4 More complete systems analysis 

using high-fidelity tools 

• Target 

Turbo-pump component analysts «♦ Entire sub-systems simulation 
Computing requirement is large; 

-♦The goal is to achieve 1000 times speed up over what was possible In 1992 


Objectives 




• To enhance incompressible flow simulation capability for developing 
aerospace vehicle components, especially, unsteady flow phenomena 
associated with high speed turbo pump. 
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INS30 - Incompressible N-S Solver 
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** Parallel version : Based on JNS30-UP 
•MPI and AALP parallel versions 

• Structured, overset grid orientation 

• Moving grid oapobility 

• Based on method of artif icial compressibility 

• Both steady-state and time-accurate formulations 

• 3»*and 5*-orderflux difference splitting for convective terms 

• Central differencing for viscous terms 

• One- and two -equations turbulence models 

• Several linear solvers : 6MRES, 65 line-relaxation, LU-56S, 

6S point relaxation, ILU(0)..,.. 

•HISTORY 

** 1982-1987 Original version of INS3D - Kwak, Chang 
** 1988-1999 Three different versions were devaped : 
INS3D-UP / Rogers, Kiris, Kwak 
INS3D-LU / Yoon, Kwak 
INS3D-FS / Rosenfeld, Kiris, Kwak 


Time Accurate Formulation 


• Time-integration scheme 

Artif icial Compressibility Formulation 

* Introduce a pseudo-time level and artificial compressibility 
• Iterate the equations in pseudo-time for each time step until 
incompressibility condition is satisfied. 


Pressure Projection Method 

• Solve auxiliary velocity field first, then enforce 
incompressibility condition by solving a Poisson equation 
for pressure. 
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Artificial Compressibility Method 
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Time-Accurate Formulation 

Discretize the time term in momentum equations using second-order 
three-point backward-difference formula 

( y*+i 

SU + dV + dW =0; = _ r .„ 

a<f dr, d( 2M 

i. * 

Introduce a pseudo-time level and artrf icial compressibility, 

Iterate the equations in pseudo-time for each time step until 
incompressibility condition is satisfied 
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Impulsively Started Flat Plate at 90° 


Aamto&xrtCtKtt 


• Time History of Stagnation Point 
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m EXP (T*n«ia, 1971) 

TNS1D-UP 

TNS3D-FS 

FmiieEkro. Sot. {Yoihida, 1985) 
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Impulsively Started Flat Plate at 90" 
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• Time History of Sta gnat ion Point 

Artificial compressibility incorporated with Poisson solver 


INS 3D -UP / DT=a0825 
EXP. (Tanedt, L971) 

- LINE-RELAXATION, NSUB=10 

■ GMRE5-1LU(DX NSUB=L0 

■ L-REL., POISSON SOL, N5UIN2 


Current Challenges 




• Challenges where improvements are needed 

- Time-integration scheme, convergence 

- Moving grid system, zonal connectivity 

- Parallel coding and scalability 

• As the computing resources changed to parallel and distributed 
platforms, computer science aspects become important. 

- Scalability (algorithmic <S implementation) 

- Portability, transparent coding, etc 

• Computing resources 

- "Grid" computing will provide new computing resources for 
problem solving environment 

- High-fidelity flow analysis is likely to be performed using "super 
node" which is largely based on parallel architecture 
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SCRIPTING CAPABILITY FOR GRID GENERATION 

> Require expertise to build scripts the first time 

> Allow rapid re-run of entire grid generation process 

> Easy to do grid refinement and parameter studies 

> Easy to try different grlddfng strategies 

> Documentation of gridding procedure 

> Written In Tel scripting language 

> works on UNIX, LINUX and WINDOWS 

> Integer and floating point arithmetic capability 

> modular procedure calls 

> easy to add GUI later If needed 
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Scripting Capability 


INPUT AND OUTPUT 


Current example: one script for each component 
(IGV, Impeller end Diffuser) 


> profile curve for hub and shroud In PLOT3D format 
(rotated by script to form surface of revolution) 

> blade and tip surfaces In PLOT3D format 

> Parameters that Can be changed 

- global surface grid spacing (on smooth part of geometry) 

- local surface grid spacing (leadlng/tralllng edges, etc.) 

- normal wall grid spacing (viscous, wall function) 

- marching distance 

- grid stretching ratio 

- number of blades 


Output 

21 > overset surface and volume grids for hub, shroud, blades 


ImbautiCttto* 


Scripting Capability 


INLET GUIDE VANES AND DIFFUSER 


Old IGV New IGV Old DIFF New DIFF 
No. of points (million) 7.1 1.1 8.0 1.6 

Time to build 1/2 day 10 sec. 1/2 day 0 sec. 

Script timings on new grids based on SGI R12k 300MHz processor 
Time to build script = 1 day for IGV, 1 day for DIFF 
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SSME-rigl /Initial start 


TIME STEP 18 / Impeller rotated 8 -degrees at 100% of design speed 
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PRESSURE 


VELOCITY MAGNITUDE 


sA SSME-rigl / Initial start 

TIME STEP 96 / Impeller rotated 42 -degrees 


vb-ooty magnitude 


- 34.3 Million Points 

- 800 physioal time steps in one 
notation. 

*One physioal time-step requires 
less then 20 minutes vail time 
with 80 CPUs on Origin 2000. 
One complete rotation requires 
one-week wall time with 80 CPUs. 
^Currently I/O is through one 
processor. Timing will be improved 
with parallel I/O since time- 
accurate computations are I/O 
intensive. With f uther 
improvements several impeller 
rotations can be completed in one 
week. 
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